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We show that isolated Weyl nodes can arise in a system of parallel topological insulator nano- 
wires arranged in a honeycomb fashion. This introduces another theoretical example of a topological 
semi-metal phase with more than one pair of Weyl nodes and due to the simple form of its Hamilto- 
nian it can be used to study various interesting phenomena associated with this phase. Our result 
emphasizes that depending on the separation of the Weyl nodes a topological electromagnetic re- 
sponse might or might not emerge and two pairs might have overall cancelling contributions to the 
net anomalous Hall conductivity. 



I. INTRODUCTION 

New topological phases of matter have been discov- 
ered recently highlighting the fact that the topological 
structures formed by the underlying electronic states are 
as important 1 as the symmetries involved and can be 
closely related to some interesting physically observable 
phenomensPtES. A recently discovered example in two 
and three dimensions is the time-reversal invariant topo- 
logical band insulator phase which can be formed by 
the electronic states in crystals with strong spin-orbit 
couplings and has been fully class ified in terms of well- 
defined topological invariant J 13 | 14 l In a strong topologi- 
cal insulator (TI) phase, the quantum states have been 
formed in a non-trivial fashion so that the system belongs 
to a different topological class than the ordinary insulator 
(01) meaning that there is no fully gapped deformation 
path in the phase space of the Hamiltonians that connects 
the Hamiltonian of the system to that of an ordinary in- 
sulator preserving the time-reversal symmetry 15 . Any 
phase space path connecting two topologically distinct 
points in the phase space would contain a point where 
the associated Hamiltonian is gapless. This intuitively 
explains the presence of the robust topological surface 
states at the interfaces between these systems and ordi- 
nary insulators. It can be understood from the fact that 
one can establish a mapping between any path connect- 
ing a point in the bulk of a topological insulator to a point 
in the vacuum and the phase space paths which respect 
time-reversal symmetry connecting two topologically dis- 
tinct regions and therefore it must contain a critical point 
where the system is locally gapless and all such points 
span the metallic surface. Therefore, any local perturba- 
tion on the surface that respects the time-reversal sym- 
metry cannot destroy the gapless surface modes. These 
gapless surface modes are then protected by the coexis- 
tence of time-reversal symmetry on the surface and the 
non-trivial topological structure in the bulk of the sys- 
tem. 

The topological protection is not restricted to just one 
and two dimensional metallic systems. Indeed one can 
generalize this idea to find a three dimensional metallic 
system where the band crossing is protected by the un- 
derlying nontrivial topological structure of the quantum 



states rather than discrete symmetries such as the inver- 
sion or the time-reversal. Recently, experimentally feasi- 
ble proposals to realize such 3D metallic phases has been 
made where energy crossings (Dirac Points) exist in the 
band dispersion and are protected by global topological 
properties of the electronic states! 16 " 20 1 In this case there 
is no discrete symmetry involved in protecting the Dirac 
point from the gap opening as opposed to the 2D case 
where we still need to keep the time-re versa l symmetry 
in order to maintain the band crossin g 21 * 22 !. The exis- 
tence of this nontrivial topological structure associated 
with the isolated energy crossings would lead to some in- 
teresting transport phenomena in the bulk and on the 
surface of such systems. A topological axion term with a 
non-uniform axion field emerges in the effective electro- 
magnetic Lagrangian which explains the chiral anomaly 
in these systems^^Sl 

In general, accidental energy crossings at certain points 
in the Brillouin zone (BZ) of a three dimensional elec- 
tronic system can only exist by fine tuning and imposing 
certain symmetries in a model. These are often hard to 
achieve and symmetries are sometimes broken in actual 
materials due to the imperfections or interactions. How- 
ever, recently it has been suggested that in certain sys- 
tems energy crossings can have topological protection. 
When such topological protection exists, the system is 
in a topological semi-metal phase even in the presence 
of weak disorder or other types of interactions that pre- 
serve the momentum conservation. These band crossing 
points in the BZ can only exist in pairs^ and are called 
Weyl points. Near an isolated Weyl point the low lying 
states can be described by a 2-by-2 Dirac Hamiltonian. 
Therefore, any local perturbation that does not violate 
the momentum conservation can only shift the position 
of such a point in the BZ and is not sufficient to gap 
out the spectrum unless it is strong enough to merge two 
such points with opposite chirality and make the system 
unstable towards becoming an insulator. These points 
can be thought of as topological defects in the fibre bun- 
dle formed by the electronic states in the BZ as the base 
manifold and are the magnetic monopoles of the pseudo- 
magnetic field associated with the gauge field defined for 
the Bloch states. Only the pairwise annihilation of such 
points with opposite chirality is possible as one can de- 
fine a well-defined charge for these monopoles which is a 
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conserved quantity in the regime where the crystal mo- 
mentum is still a good quantum number. 

The question is how such energy crossings can occur in 
a band system. Theoretically, it is possible to choose a 
physical parameter to adiabatically drive a system from 
a topological insulator phase to an ordinary insulator or 
vice versa by tuning such a parameter. As one fine tunes 
such a parameter to the critical point where the gap 
closes (which is inevitable when we have time-reversal 
symmetry) the system becomes an unstable bulk metal if 
the time-reversal and inversion symmetry is preserved at 
the same time. At this critical point the system would be 
in an unstable 3D metallic phase with degenerate bands 
near the crossing points. This gapless metallic phase is 
susceptible to local perturbations and instabilities that 
drive the system to either topological or ordinary insula- 
tor phases. Now if another parameter in the Hamiltonian 
can be tuned in such a way that it separates the degener- 
ate bands at the crossing points in the BZ (by breaking 
the inversion and/or the time-reversal symmetry) then 
the system would be in a topological semi-metal phase 
with gapless modes dispersing in three spatial dimen- 
sions which are "robust" against momentum conserving 
perturbations. The topological nature of the phase re- 
flects itself in the appearance of a nontrivial term in the 
effective Lagrangian of the electromagnetic fields^. This 
topological term is an abelian axion term with a non- 
homogeneous axion field 



Saxion = C J 6(X)TAT, (1) 

where C is a constant and 6{x) modulates in space and 
time with wave- vectors that depends on the separation of 
the Weyl nodes in momentum and energy and for a single 
pair of Weyl nodes it is given by 0(x) = q-x in which q is a 
four- vector of the separation of the pair of the Weyl nodes 
in the crystal momentum and energy^"^. This topolog- 
ical term in the effective electromagnetic Lagrangian is 
closely related to the topological transport phenomena 
that is present in these systems^. It turns out that a 
system with more than one pairs of Weyl fermions might 
have zero net anomalous Hall conductivity since pairs can 
have total cancelling contributions due to the exactly op- 
posite separations of the them in momentum space. Here 
we consider a simple theoretical model which can lead to 
a topological semi-metal phase with two Weyl pairs under 
certain circumstances. Although there are big challenges 
on the way of the experimental realization of such a sys- 
tem, it can be used as a platform to study the bulk and 
the surface properties of a system in a Weyl semi-metal 
phase when there are more than one pair of Weyl nodes. 

The paper is organized as follows: In section II, we 
provide a simple mathematical proof of the topological 
protection of the Weyl nodes using the Clifford algebra 
when we have a pair of them separated in the BZ. In sec- 
tion III, we introduce a lattice model made by arranging 
parallel TI nano wires in a honeycomb fashion with a lat- 



tice spacing which is small enough to allow a significant 
hopping of the electrons between these wires and then 
we discuss how one can achieve a topological semimetal 
phase out of this arrangement. 

II. TOPOLOGICAL SEMI-METAL 

At low energies, the minimal Hamiltonian describing 
the unstable 3D gapless system at the TI-OI transition 
point with a pair of energy crossing points would have 
the following form 

n =J2 E *£(*)[Aa.(fc-tfl„)]* a (fe), (2) 

a=± \k—w a \<A 

where X a are 2-by-2 matrices which are given in terms 
of Pauli matrices, cr$, by 

A = (v x a x , v y Gy, v z a z ), (3) 

w+ = w = wq are the position vectors of the Weyl 
points with +1 and —1 chirality ( = sgn(v x v y v z )). This 
Hamiltonian consists of doubly degenerate Dirac cones 
with the Dirac points located at wo. The two compo- 
nent spinor operators ty± are associated with the modes 
on the Dirac cone with chirality a = ±. This Hamilto- 
nian describes the low-energy gapless phase at the phase 
transition point between a topological insulator and an 
ordinary insulator. The cutoff, A, defines the momen- 
tum range around each Weyl point for which the lin- 
ear approximation is valid. When the opposite chirality 
Weyl points happen at the same point in the BZ, there 
is no topological protection and even momentum con- 
serving perturbations that connect the degenerate states 
can drive the system to an insulator phase. Now, if 
somehow the opposite chirality points become separated 
{w + —W- = q 7^ 0) then the system would be in a gapless 
topological phase where small local momentum conserv- 
ing perturbations cannot produce a gap in the spectrum 
and would only shift the position of the Weyl points in 
the BZ. 

Using the mathematical properties of the Clifford alge- 
bra, it is possible to see how such protection can occur for 
a gapless system containing a pair of separated opposite 
chirality Weyl points. The local Hamiltonian matrix in 
the BZ describing the low energy states in the topological 
semi-metal phase containing a pair of opposite chirality 
Weyl points can be written in terms of 4-by-4 Dirac ma- 
trices as follows 

H{k) = k xll + k yl2 + k z <y 3 + |t£, (4) 

where for simplicity and without the loss of generality 
we have put \v X}VyZ \ = 1. The first three 4-by-4 matrices 
7i,2,3 are mutually anti-commuting to ensure that we will 
get two opposite chirality Dirac nodes dispersing in three 
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dimensions. Without 73, we get two degenerate Dirac 
cones, therefore we must introduce it into the Hamilto- 
nian to separate the Weyl nodes by q in momentum space 
and it must commute with at least one of the other ma- 
trices or a momentum independent linear combination of 
them depending on the direction of the desired separa- 
tion. This commutation relation ensures that the effect of 
the added term is only opposite shifts of the Weyl nodes 
along an specific direction in the momentum space. We 
consider Weyl nodes separated in k z direction, therefore, 
we must choose 7 3 in a way that it commutes with the 
matrix associated with k z , i.e., 73 and anti-commutes 
with the rest. Note that by writing the Hamiltonian 
local in momentum space, we have implicitly excluded 
any perturbation that connects two Weyl nodes which 
might arise from the scattering sources. Those can po- 
tentially gap out the Weyl nodes by connecting states 
with different momentums, however, the irrational sepa- 
ration of the Weyl nodes ensures that the system is not 
unstable toward the charge density wave ordering with 
a commensurate wave vector. To summarize, we have 

{liilj} = {71,2,73} = and [73,73] = °- Note that 
the low energy Hamiltonian describing any pair of oppo- 
site chirality Weyl nodes is connected by a well defined 
unitary transformation to the Hamiltonian given by Eq. 
Q. The spectrum of this Hamiltonian can be obtained 
by, first, multiplying the Hamiltonian by itself and then 
taking the square root of the resulting diagonal matrix. 
It is easy to check that the relation between the above 
matrices leads to the following spectrum for this Hamil- 
tonian 



e Ks = ±^kl + kl + (k z + ff, s = ±l, (5) 

and it has all the correct properties, i.e., the opposite 
chirality Weyl nodes separated by q in the k z direction. If 
there exist a momentum conserving perturbation which 
opens up a gap in the spectrum of the above Hamilto- 
nian, there must be an associated hermitian matrix that 
anti-commutes with all of the terms in T~L(k). If such a 
term does not anti-commute with all the terms, then it 
can only change the position of the Weyl nodes. This 
can be seen by following the same steps we have taken 
to obtain the spectrum of H(k). The statement we want 
to prove here is that no such term exists and any 4-by- 
4 Hermitian matrix which can be written as a super- 
position of the normalized Clifford algebra matrices, 7^, 
(j g 2 = ^4x4) commutes with at least one of the terms in 
this Hamiltonian. Therefore, it cannot produce a gap in 
the spectrum since this way it can only shift the posi- 
tions of the Weyl nodes in momentum space. In order to 
prove this statement, we show that the maximum num- 
ber of mutually anti-commuting matrices for which one 
of them commutes with 73 is three. In order to do so 
we need to choose a representation of the present ma- 
trices. Since 73 and 73 commute with each other, we 
can choose a basis in which both of these matrices are 



diagonal. This way we have 73 = diag[l, ai, a 2 , — 1 — 
ai - a 2 ] and 7 3 = diag[l, 61, 6 2 , -1 - bi - b 2 ). It 
turns out that only three independent choices can sat- 
isfy this commutation condition, i.e., (ai, a 2 , 61, 6 2 ) = 
(1, -1, -1, 1), (-1, 1, 1, -1), (-1, 1, -1, -1). For each of 
these choices we can show that only two independent 
normalized matrices exist which can simultaneously anti- 
commute with 73 and 73 and with themselves. This 
would then imply that the maximum possible number 
of mutually anti-commuting matrices, when we require 
that one of them to commute with 73, is three . These 
three have already been exploited in the Hamiltonian by 
coupling to three momentum components and therefore 
no other matrices can be added to anti-commute with all 
of these matrices and 7 3 in order to gap out the spectrum. 
To see how it works, consider the most general normal- 
ized matrix satisfying the aforementioned relations for 
the first choice of (ai, a 2 , 61, 62)- For this particular rep- 
resentation it has the following form 
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(6) 



Matrices of this form that anti-commute with each 
other must satisfy |^ — \ = \6i — Oj \ = tt/2. It is easy 
to see that only two independent matrices can be found 
in this representation to satisfy these equations. These 
matrices are already exploited in the Hamiltonian by cou- 
pling to momentum components k x and k y and therefore 
no other matrices can be added to the Hamiltonian given 
in Eq. Q that anti-commutes with all the present matri- 
ces. The same argument can be established for the other 
two choices of (ai, a 2 , bi, 6 2 ). Therefore, Weyl nodes can- 
not be annihilated unless two of them with opposite chi- 
rality merge at the same point in the BZ. This cannot 
be possible without the sufficiently strong local pertur- 
bations (~ \hvpq\ to bring two nodes to the same point 
in BZ) and the system is in a topological semimetal phase 
in at least a region of the phase space. The above proof 
is a generalization of the previously known argument in 
the literature which one considers only one of the isolated 
Weyl nodes in a two band model using two-by-two Pauli 
matrices. Here we have considered all the four bands 
which are present in the low energy theory describing a 
pair of Weyl fermions. In the rest of this paper we present 
a system with more than one pair of Weyl nodes where 
the stability arguments should be considered more care- 
fully as it might be possible to connect two pairs of Weyl 
nodes with a charge density wave ordering. 



III. TI NANOWIRES 

The surface modes of a cylindrical topological insulator 
system can be described by a Dirac Hamiltonian in a 
curved space^ (h = 1) 




FIG. 1: (Color online) The surface spectrum of a wire made 
from a strong topological insulator in the presence of half 
integer multiple of the magnetic flux quantum along the wire. 
The bandwidth in which the spectrum is nondegenerate is 
A = hvo/R. Where R is the radius of the wire. The spin 
texture is shown schematically for the states dispersing along 
k z with the blue arrows. 
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FIG. 2: (Color online) Two adjacent parallel wires. The elec- 
tronic states on one can interact with those on the other. 



when there is a magnetic flux through the wire equal to 
the odd multiple of the half quantum of the magnetic 
flux r] = n + 0.5 we get a non-degenerate gapless band 
for I = n 



H k = ^ I V • ™ + h ■ (P x °) + (P x °) ' n]> ( 7 ) 

where a is the pauli matrix vector which acts on the spin 
Hilbert space, p = — zV is the momentum operator and 
n is the unit vector normal to the surface. For a cylinder 
of radius R along z axis we have ft = cosipx + sincpy. 
The Hamiltonian that governs the TI surface modes in 
the presence of a sufficiently thin magnetic flux, <j> = 77^0, 
can then be written as 



We can simplify this Hamiltonian by a unitary trans- 
formation Hk = U^(cp)HkU(cp) in which U(tp) is given 

by 

using this transformation we can get rid of the phase in 
the off-diagonal components of the Hamiltonian matrix 

Hk ~[ ^ _i( i ^ + l/ _|) > ). 

The eigenstates of this Hamiltonian are given by 
i/jklip) = ^ l<pl ^ki where ipki are eigenstates of the H^i 
defined as 

H k i = k z a 2 - ^ a s(l + ^ - rj), (11) 



= k z * 2 , (12) 

We assume that R is small enough in a way that there 
exist a significant energy range, A, in which the gapless 
band does not overlap with other bands (See Fig. [I]) then 
we can use this energy band as a building block of the 
model to realize a Weyl semimetal phase. 



IV. LATTICE OF PARALLEL WIRES 

When the distance between two parallel nano wires is 
sufficiently small, there would be an overlap between elec- 
tronic wave-functions and this leads to a hopping between 
adjacent electronic states of the surface modes. Here we 
consider a honeycomb arrangement of these nano wires 
by considering the nondegenerate gapless modes of these 
wires that can be achieved by exposing them to a static 
magnetic field along the wires. We assume that the spin- 
momentum locking that happen at the surface of topo- 
logical insulators are in opposite directions for wires in 
the A and B sublattices. We should point out here that 
all the by far discovered topological insulators happen to 
have the same direction for the spin-momentum locking. 
In theory, one can get a system with an opposite spin- 
momentum locking by simply changing the sign of the 
spin-orbit coupling required to get the original topolog- 
ical insulator system^. We also assume that the Dirac 
points for both A and B wires have the same energy, how- 
ever, in the end we relax this assumption and we discuss 
how this is essential in order to get a Weyl semimetal 
phase. On the other hand, these wires are multi-band 
systems, therefore, they must be sufficiently thin to allow 
a significant gap, A, in which the gapless band does not 
overlap with other subbands (see Fig. With these 

assumptions and by considering only the non-degenerate 
lowest lying band in each wire in the energy range 2A 
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(see Fig. [T]) we can split the Hamiltonian into two parts 
as follows 



H = H + H t 



(13) 



in which Ho which sums up the gapless nondegenarate 
modes of each individual wire in a second quantized no- 
tation and it is given in terms of the Fourier components 

as 



J *E*liW(^)* fcl w, (14) 



where T3 is the pauli matrix acting on the sublat- 
tice Hilbert space, k z is the momentum along the wire 
(Fig. [T]) and k± spans the honeycomb's reciprocal lat- 
tice (similar to what one gets for graphene) and Sfr^ = 

C0a+> ^A-' ^B-)- ^ ne fi rs ^( ^ as ^) two components 
act on the spin space of the wire in the A( B) sub-lattice. 

^A/Ba crea ^ es an electron in the a branch of the A/B 
wire when it acts on the vacuum state. Using this rep- 
resentation, H t which represents the direct hopping be- 
tween nearest-neighbour wires, can be written as 



Hi 



' J #„ 











>1 



in which g(kj 



-tJ2 5 - e~ ik ^ Si and S 2 



(15) 

are three vec- 



tors that connect a site in honeycomb lattice to the three 
adjacent points. The energy spectrum of the total Hamil- 
tonian would then be doubly degenerate 



e±(k z ,k±) = ± 



vlk 2 z 



(16) 



g(k±) vanishes at two inequivalent points in the two- 
dimensional reciprocal lattice, i.e., at and is lin- 
ear near these points. Therefore, we get two degenerate 
three-dimensional Dirac points at (k z , k±) = (0, K^). In 
order to get a topologically protected semimetal phase 
we need to break time-reversal or inversion symmetry in 
such a way that it separates the Dirac points. Breaking 
inversion symmety can be easily realized by relaxing the 
assumption we had in the beginning, i.e., Dirac point at 
two sub-lattices now can have different energies. By re- 
vising Ho to account for such an energy difference, V, we 
get 



Ho= I ^^^(^)(^r 3 a 2 + W 3 )^(^), 

(17) 

In this case the degeneracy is lifted and we have four 
bands £ sr (k Zl k±) (s, r = ±) given by 



-(**, k ± ) = ry/(v k z + sV)2 + \g(k ± )F 



(18) 



K+ 




K 



±A ■ k 




FIG. 3: (Color online) The three dimensional BZ before and 
after separating the Weyl nodes along the k z . We have two 
inequivalent pairs of opposite chirality Weyl nodes. 



This spectrum has two pairs of Weyl points in the three 
dimensional BZ at which the bands cross linearly. Each 
pair consists of two opposite chirality Weyl points cen- 
tred at (k z ,k±) — (0,-R^). They are separated along 
k z axis by Q = 2V/vo- According to the argument pre- 
sented in the previous section, these crossings are topo- 
logically protected against various local momentum con- 
serving perturbations as long as they are separated in 
the BZ. It is important to note that although these iso- 
lated Weyl nodes are robust against local momentum- 
conserving perturbations, the peculiarity in the Hall re- 
sponse that is present in a system with only one pair of 
Weyl nodes does not exist here, i.e., the transverse con- 
ductivity in the presence of the applied magnetic field 
along the wires, a xy is zero 1 16 * 31 ! and two pairs have can- 
celling contributions to the Hall conductivity as they are 
separated in the opposite way at two valleys consider- 
ing their chirality. There is also a possibility of phase 
transition to an insulator phase due to the electron- 
electron interactions. This instability arises in a mean- 
field treatment of the inter-rod electron-electron interac- 
tions. When the wire radius, R, is much smaller than the 
honeycomb lattice constant, a, the coulomb interaction 
between electrons in the adjacent wires can be written as 



H e 



U 



/L/2 nC 
dz 
-L/2 J-i 



du 



h(R, z)n(R + 5,z + u) 



(19) 

where n is the electronic density operator and U = 
e 2 /(47rea). R, S and z have been defined in the Fig. 
pi In a mean-field treatment of the above term and by 
considering a Kekule type modulating order parameter 
with a wave- vector G = K + — K~ , it is possible to con- 
nect the Weyl points separated by G. This would then 
open up a gap and the system becomes an insulator. The 
critical coupling, U c at which such a phase transition to 
an insulator phase occurs can be computed for this sys- 
tem. It is a function of the potential difference, V, as 
well as the hoping strength, t, and is given by 
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which is similar to the term which has been introduced in 



FIG. 4: (Color online) The modulations in the magnetic flux 
with a zero average through each hexagon introduces a com- 
plex next nearest neighbor hopping parameter, t '. The 
vectors that connect adjacent sites in the same sublattice are 
±(&i, 62, 63) = ±(<5i — £2, 82 — £3, £3 — Si) 



U c =- h+m— 



(20) 



in which a = (8/V^tt) 1 ^ 2 . This critical coupling would 
be of the order of the single wire's non-degenerate band- 
width A (see Fig. [TJ whenever the hoping strength is 
significant enough. For U < U c the system would re- 
main in the topological semimetal phase and for U > U c 
the system would become a three-dimensional version of 
Kekule insulator which has been discussed previously in 
two dimensions for graphene^. It is possible to make 
the system stable against such phase transition for all 
ranges of interactions at least in this mean-field channel 
by introducing a term in the Hamiltonian which sepa- 
rates two Weyl nodes at each valley (K + and K~) by a 
different amount. In this case the system becomes stable 
against the perturbation that connects two valleys since 
all the opposite chirality Weyl pairs would have incom- 
mensurate separation in momentum space after addition 
of such a term. This term can be induced by consid- 
ering next nearest-neighbour inter-wire hopping of the 
electronic states. The hopping amplitude is imaginary 
considering a modulating magnetic flux with a zero net 
average through each hexagon (see Fig. 4). This can be 
achieved by applying a modulating magnetic field instead 
of the uniform magnetic field that is required to induce 
half-quantum magnetic fluxes in the wires. Such a term 
was first introduced for a honeycomb lattice to realize a 
quantum Hall phase in a system without a net uniform 
applied magnetic field^. The Hamiltonian for this next- 
nearest neighbor hopping in the momentum space can be 
written as 



H' 



k± 

(21) 



Eq. (17) but here the V which separates the Dirac nodes 
along k z is a function of the k± and takes different values 
at K + and K~ . The presence of /x(fej_) separates two 
pairs in energy since it also depends on k± and can take 
different values at two valleys. In the original Haldane's^ 
flux configuration /i(fcj_) and V(k±) are given by 



V(k±) = — 2£'sin0 ^ sink_\_-bi 



2=1,2,3 



and 



fjb(k±) = 2t f cos <p cos k± ■ bi 



(22) 



(23) 



=1,2,3 



bi and <p have been defined in the caption of Fig. 4 and 
t' is the next-nearest neighbor hopping amplitude. 

Therefore, the distance between each pair of Weyl 
nodes is now an incommensurate wave-vector since the 
separations of Weyl nodes at two valleys, i.e., qi = 
2V(K + )z/vq and qr 2 = 2V(K~)z/vq are not exactly 
opposite. This makes the system stable against charge 
density perturbations that connect quantum states of two 
valleys. The Axion field is not zero in this case as the 
separations are not exactly opposite and instead it mod- 
ulates in spaced 

Therefore in the presence of the Haldane term the sys- 
tem breaks the time reversal symmetry and it would have 
a nontrivial topological response. In this case the trans- 
verse conductivity in the presence of the magnetic field in 
the z direction can be obtained by a summation over the 
transverse conductivity of the m assive two-dimensional 
Dirac fermions for each fcj 16 * 31 ! . At k± = K + and 
—Qi<k z <qi the sign of the masses of two-dimensional 
Dirac states are in such a way that they contribute con- 
structively to the transverse Hall conductivity. Similarly, 
at k± = K~ and —q 2 < k z < q2 the transverse Hall con- 
ductivity is nonzero. The net o xy can be obtained, when 
the chemical potential is inside the gap in the presence 
of an applied magnetic field in the z direction, by a sum- 
mation over the contribution of all two-dimensional Dirac 
fermions labeled by k z and the valley index, i.e., 1 and 2 



cr xy = / Wi(k z ) - i/ 2 (k z )] 

J-A 



dk z e 2 



(24) 



in which v\ = &(qi — \k z \) and 6(^2 — \k z \)- Using Eq. 
( 22 ) the transverse conductivity becomes 



6^' 



sin ( 



nh 



(25) 



This result highlights the fact that although the mo- 
mentum conserving perturbations cannot gap out the 
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system, this does not necessarily imply the existence of 
a topological electromagnetic response in the system and 
two pairs of Weyl nodes can have cancelling contribu- 
tion to the electromagnetic response in special cases and 
therefore (7 xy can be zero even when the Weyl nodes 
are well separated in the BZ which is the case for our 
model in the absence of the Haldane term. Adding the 
next nearest hoping and introducing the Haldane term 
would lead to a nonzero net topological electromagnetic 
response since two Weyl pairs would now have different 
separations and therefore their contributions to the trans- 
verse conductivity do not sum up to zero. 

V. CONCLUSIONS 

To summarize we have considered a system of parallel 
nano wires in the time-reversal invariant topological insu- 
lator phase and by breaking inversion and time-reversal 
symmetry we found that it is possible to realize a topolog- 
ical semimetal phase with two pairs of Weyl nodes in the 
three dimensional BZ which are protected against local 
perturbations that conserve crystal momentum. A non- 
trivial topological electromagnetic response might arise 
under certain circumstances when the contributions from 
two pairs does not exactly cancel each other. It would 
also be interesting to study the surface modes in the pres- 
ence of more than one pairs of Weyl nodes which can be 
done using the model introduced in this paper. 

Finally, we point out that although the growths of 
large-scale vertically aligned nano-rods and nanopillars 
with various arrangements have already been achieved 
by experimentalists^^, realizing a system for which one 



can observe the robust Weyl energy crossings seems not 
to be experimentally plausible at this point. One can 
name various obstacles on the way of its experimental 
realization. First, wires are multi-band systems and the 
two-band model approximation used throughout this pa- 
per requires a significant gap for the other existing sub- 
bands in the nano- wires. For the so far discovered topo- 
logical insulators, the surface Fermi velocity^, vq, is of the 
order of 10 5 m/s (5 x 10 5 m/s for Bi 2 Sefl), therefore, the 
wire diameter required (2R < 60nm) to get a significant 
gap (A > lOmeV) is still out of the experimentally feasi- 
ble sub-micron ranges^"^. Finally, another challenge to 
realize such a system is to find two types of topological 
insulators with opposite spin-momentum lockings which 
although in theory nothing prevents having such systems, 
in all the so far discovered topological insulators surface 
electrons' spins curl around the Dirac point in the mo- 
mentum space in the same directions. We note that the 
direction of the spin-momentum locking is a material de- 
pendent property and in the theoretical lattice models 
for the topological insulators in three dimensions, it is 
possible to change it on the surface by varying physical 
parameters of the bullP^] 
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